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Abfltntct 


The isotopic compoeitiona of galactic cosmic ray boron, carbon, 
and nitrogen have been, measured at energies near 300 MeV amu~^ using 
a baUoon-borne instrument at an atmospheric depth of ~ 5 g cm~^. The 
calibrations of the detectors comprising the instrument are described. 
The saturation properties of the cesium iodide scintillators used for 
measurement of particle energy are studied in the context of analyzing 
the data for mass. The achieved rms mass resolution varies from 
~0.3 amu at boron to ~0.5 amu at nitrogen, consistent with a 
theoretical analysis of the contributing factors. Corrected for detector 
interactions and the effects of the residual atmosphere, the results are 
*08/8 = 0.3318:11, *0C/C =0.06l8;i?, and *0N/N = 0.4218;]?. A model of 
galactic propagation and solar modulation is described. Assuming a 
cosmic ray source composition of solar~like isotopic abundances, the 
model predicts abundances near earth consistent with the 
measurements. 
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Chapter 1 
tatroduotton 

Within the last decade, increasing* emphasis in cosmic ray 
astrophysics has been placed on the measurement of the isotopic 
composition of nuclear species with atomic number Z2 3, for several 
reasons (e.g. Stone 1973). Because the cosmic rays are extremely young 
(~ lO’’^ yr) in comparison with the solar system and constitute the only 
material from outside which is directly observable, it is of interest to 
know whether the astrophysical conditions under which they were 
synthesized is similar to that of solar system material. Measurements of 
the cosmic ray elemental composition (e.g. Garcia-Munoz and Simpson 
1979, Lezniak and Webber 1976), although experimentally more 
tractable, do not always bear directly on the question of source 
composition because of possible Z-dependent selection effects (Casse' 
and Goret 1978) on the acceleration of the particles subsequent to their 
synthesis and injection into the interstellar medium. Also, because the 
cosmic ray abundances observed near earth are contaminated by the 
products of nuclear interactions suffered in the interstellar medium en 
route to the solar system, and the cross sections for such interactions 
are mass~dependent, only by measuring the isotopic composition can 
one adequately deconvolve observed abundances into source abundances 
and at the same time understand the nature of the galactic propagation 
process. Finally, precise measurements of the cosmic ray age are 
possible if radioactive nuclides with half-lives of the order of that age 
(e.g. ^^5e, ^^Al) can be resolved from the more abundant neighboring 
isotopes. 
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Th« elcmenU boron, carbon, and nitrogen, and the ieotopei of 
each, have different hietoHes regarding their mode of origin. Boron hae 
onlj a transient existence in iitars as it is extremely unstable at 
temperatures required for the synthesis of elements heavier than 
helium. That its abundance reUtive to carbon in the cosmic rays (as 
well as that of lithium and beryllium) is orders of magnitude above the 
solar system value, is evidence that cosmic rays trd.verse several g cm~^ 
of interstellar matter before they are observed near earth. A 
measurement of the quantity ^°B/B in the cosmic rays depends to 
ffrst-order only on the ratios of the relevant cross sections for 
producing *®B and **B from the spallation of heavier species, primarily 
carbon and oxygen. 

The moat abundant isotope of carbon (^^C), on the other hand, is 
copiously produced in stars via helium burning. The rarer isotope is 
produced as one of the products of the CNO bi-*cycle in hydrogen 
burning. Its abundance, when the cycle is operating in equilibrium, is a 
function of the temperature of the star. The solar system value of 
i^C/C w 0.011 contrasts somewhat with radio observations of molecular 
clouds where *®C/C is a factor of ~1.5 higher (Wannier 1980). At 
temperatures > 10^ one can expect a value as high as ^^C/C w 0.2, 
which, in fact, has been observed in some carbon stars. It should be 
noted, however, that using cosmic ray measurements of *®C/C to infer a 
source abundance is difficult because a large portion of the observed 
is secondary (results from the spallation of N, 0, etc.). 

The production in stars of the isotope is also associated with 
the CNO bi-cycle. In fact, if the cycle reaches equilibrium, essentially all 
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of ihtt s««d C and 0 nuolai are oonrerted into The iaoiope ii al*<o 
produced, but ita abundance if predicted to be **N/N 4x 10~®, 
approximatelj two ordere of magnitude below what if obeenred in the 
flolar f3rftem. Larger amounta of might be produced in ezploaive 
CNO burning (Truran 1977). Lihe carbon, inferring a aource abimdance 
from coamic ray measurementa of *®N/N requlr'*n reaaonably preciae 
valuea of the croaa aectiona to produce and from abundant 
heavier t;peciea (predominantly ^^0) ao aa to accurately account for the 
aecondary component of the obaerved flux. 

Recent coamic ray meaaurementa of *®C/C, and *®N/N can 

be divided into two categoriea, eaaentiaUy on the baaii of energy. In the 
first category are relatively !on energy (~ 100 MeV amu”*) satellite 
measurements (e.g. Garcia-Munox sf aL 1977a, Guzig 1980, Wiedenbeck 
9t al, 1979). Although the satellite instruments typically have 
reasonable-to-excellent mass resolution capabilities, the interpretation 
of the measurements is complicated by the fact that the relevant cross 
sections at the corresponding energies in interstellar space vary 
considerably with energy, and measurements of some of the more 
important ones (e.g. ^*0+p-+ *®N,^*N) do not exist. In the second category 
are balloon measurements at energies of '^300-400 MeV amu"* 
(e.g. Hagen et al 1977, Buffington vt al. 1978, Webber and Kish 1979, 
Webber af. 1979). Disagreements among the results of these 
experiments, which are apparently statistically significant, suggest 
instead that systematic effects in the measurements may not have been 
properly accounted for. An advantage of the higher energy 
measurements is that the cross sections vary less with energy, emd that 
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th« uie of ■emi'^mpirioal formulao (Silbarberg and Tiao 1973a, 1973b, 
1977a, 1977b), which are baaad lu , >irt on high energy measurement! 
(Lindstrom wl * > 1976) is likely to be more reliable. 

We present in this thesis the results of new measurements, using a 
balloon'-bome detector array of position-’sensitive proportional 
counters and energy^mensitive cesium iodide (Csl) scintillators, of the 
isotopic compositions of the elements boron, carbon, and nitrogen in the 
cosmic rays, in Chapter 2 we describe the technique and the Instrument 
used in making the measurements, and discuss quantitatively the 
calibrations of the detectors. We go on in Chapter 3 to show in detail the 
steps involved in calculating masses, and in doing so deduce a 
light-energy relationship for Csl. The resulting mass distributions are 
analyzed for fractional abundances using a maximum-likelihood 
technique. We present a detailed analysis of the factors which 
contribute to the mass resolution and compare the calculation with the 
performance of the instrument. In Chapter 4 we correct the 
measurements for the contamination due to the residual atmosphere 
above the detector. We then describe a model of galactic propagation 
and solar modulation which allows one to predict abundances observed 
at earth given a set of source abundances. Finally, we compare our 
results and those of other recent measurements with the predictions of 
the model. 


\ 
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Chuptinr 2 

bntnuaentKtion ind Calibratioiui 


2,1 The Inairuxnent 

The Caltech High Energy Isotope Spectrometer Telescope (HEIST) 
is a balloon*~bome cosmic ray detector array designed to measure 
charge, mass, and energy of individual costnic ray nuclei with nuclear 
charge Zk 4. The instrument has been described in detail by Wiedenbeck 
(197B). We review here the major features. 

In Figure 2.1 we depict the basic technique employed by HEIST. A 
cosmic ray nucleus of mass M [proton mass units (pmu); 
1 pmu" 1.0073 amu], atomic number Z, and total kinetic energy E 
traverses a detector of thickness t at an angle d in which it loses an 
energy AE, and stops in a subsequent detector where it loses its 
remaining energy E'»E-~AE. Since the range R in a given material is a 
function of energy, mass, and charge we have 

R(E'+AE,M,Z) * t seci>+ R(E',M,Z) (2. 1 ). 

To understand how R varies with E, M, and 7, we consider the 
energy loss per unit pathlength (Rossi 1952, Jackson 1975, Janni 1966) 

=0.307MeVcm2g-*-?2?![in(^E£^Z!£)-^2^C-|] (2.2) 

pdx ® I ' ^ 2 -’ 

where 

s particle velocity (in units of the speed of light c), 

7 = (1 

Zm« Ani = mean atomic number, weight of detector material, 

p = density of material (gcm"'^), 

m*c* = electron rest mass, 0.51 1 MeV, 

I - mean ionization potential of material, 

C = correction for atomic shell structure, 

A 

— = correction for density eflect. 


S 
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Fifore 2.1 

The AE-E' technique of isotope identification. As the particle enters the 
AE-detector it has a total kinetic energy E. As it enters the E’-detector 
it has an energy E*. Since the pathlength in the AE-detector is tseci?, ire 
have R(E'+AE.M,Z)»R(E',M,Z)-»-tsecn>, from which we can extract M. 
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Since velocity fi depends only on the kinetic energy per unit mass E/M, 

dP o R 

we thus have for a given material « Z*S(~), where S(e) is the specific 
ionization of a proton of kinetic energy e. The range is 


R(E,M,Z) » / 


dE 


b dE/dx 


where Rp(e)« / de-^ 


S(e) 


is the range of a proton of kinetic energy e. 


Although the 


J1 

Z2 


scaling of range begins to break down for 


2Z 

137 


E 

corresponding to — ^ 4 MeV aniu"' for Z=8, it is valid for energies and 
M 

charges of interest in this experiment. 


For illustrative purposes we use a power law approximation 
Hp(c) = kc* (accurate to ±15% for Kc< 1000MeV; depending on the 

material, a is t 3 rpically 1.7) so that R(E,M,Z) = •^^(‘^) • this 

Z* M 


approximation, Equation 2.1 can be solved for M explicitly: 
M = 

1 ? t sect? 


(2.4). 


Thus knowledge of Z and t together with measurements of AE, E'-E-AE, 
ahd -d suffice to determine M. We defer a more detailed discussion of the 
technique to Chapter 3. 


Our instrument consists of two major components (Figure 2.2). 
First is a hodoscope consisting of eight multiwire proportional counters 
(MWPC's) X1-X4 and Y1-Y4 for measurement of particle trajectory. 
Anode wire spacing is 4mm, cathode wire spacing is 2mm. A mixture of 
70% Ar and 30% COg at a pressure of 1 atmosphere is used for counter 
gas. Readout of individual counters is by means of an electromagnetic 
delay line (~14nsecmm"^ delay) to which the cathode wires are 
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Fiftune 2.2 

Elements of the -HEIST detector (from Wledezzbeck 1970). 

a) Cross sectional view of the instrument. The cross section is taken 
along a diagonal of the proportional counter hodoscope. The light pipes 
and attached photomultipliers lie off of this plane and are shown 
projected onto it. 

b) View from above. 

c) View taken through section A-A indicated in b). Note that this view is 
at 45® to that shown in a). Each of the detectors DO through D7 has 
discriminators corresponding to whether the output from PMT dynode 3 
(H for high level) or dynode ? (L for low level) is used. DO and D1 each 
have in addition a "medium level" discriminator (M;. The stopping 
equation has been designed to eliminate 2 particles when possible but 
to accept nuclei with 4. 
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c 


XI 

Yl 

X2 

Y2 


X3 

Y3 

X4 

Y4 



01 (3mm) 
03 (4mm) 
09 (8mm) 
07((7mm) 


Coincidence equations: 

Stopping = (XI + X2)-(YI+Y2)-(X3 + X4)*(Y3+Y4)* 

[(DOL-DIL+ DOL -D2L + 01 L • D2L ) • (D4L+ DSL) + DOM • Dl M ] • D8M • A 

Penetrating = (XI -•■X2)-(YI+Y2) -(X3+X4)-(Y3+Y4)- 

(DOL-DIL + DOL-D2L + DIL-D2L)-(D4L+D5L)*D8M 

where L = low level, M = medium level 
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capacitively coupled. Each coordinate li divided into three sectors of 
176 mm width. Comparinona of the signal amplitudes among sectors 
determine the gross position of a particle, while timing of delay line 
pulses determines the position of a particle within a sector (fast 
electronics capable of 1 nsec resolution are used for delay line timing). 
The active area of each MWPC is ~ SOcmx 50 cm. 

For measurement of particle energy we use a stack of eight 
cylindrical CsI(Tl) scintillation crystals, DO through D7, each nominally 
11 inches (27.9 cm) in diameter, ranging in thickness from 3 mm to 
17 mm. Each cr)rstnl is sandwiched between two sheets of 3.3 mg cm"® 
millipore to minimize absorption of scintillation light at the crystal 
faces. Scintillators are optically separated from each other by 6.4 /im 
aluminized Mylar. One~balf of the circumferential edge of each of the 
crystals is optically coupled to a Lucite light pipe, by which scintillation 
light from the crystal is transported to its own 5" photomultiplier tube 
(PMT). Associated with each PMT are analog and digital circuits for pulse 
height analysis of PMT output. To maxisuze the dynamic range the 
outputs from both the third (’high level" or 'low gain') and seventh ('low 
level" or 'high gain") dynodes of each PMT are amplified and shaped 
independently, and discriminators direct one of these signals to an 
analog-to-digital converter for readout. D6, a ninth Csl crystal used to 
tag particles which penetrate D7, has only discriminators (is not pulse 
height analyzed). Finally, two NE 102 plastic scintillators are mounted 
below the Csl stack and used in an anticoincidence mode to tag wide 
angle penetrating events which enter and leave the Csl before 
encountering D8. 
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The initrument ie designed to operate in either of two modes 
called '^topping" and 'penetrating*'. The detailed coincidence equations 
are given in Figure 2.2c. Stopping events are simply those for which the 
particle stops somewhere between the top of D1 and the bottom of D7, 
and are suitable for our AE-E' mass anal 3 rsis technique described above. 
We confine ourselves in this thesis to mass anal 3 rsis of particles which 
stop in D7 (Chapter 3). Penetrating events are those for which the 
particle penetrates beyond the bottom of D7. Because we cannot 
analyze these events for mass, they are accumulated in fiight with a 
lower priority than stopping events and are useful mainly for in-flight 
calibrations. Both stopping and penetrating events require that the 
particle trigger a sufficient number of MWPC’s to adequately deduce 
particle trajectory. 

The payload was launched from Yorkton, Saskatchewan, Canada at 
0303 CST on 30 August 1978. The magnetic rigidity cutofi at this 
location is ~700MV/c, corresponding to c~65MeVpmu"^ Prior to 
flight, the electronics associated with MWPC X2 failed and were tux'ned 
ofl^. Approximately twenty hours of data were taken at an atmospheric 
depth of 4-6 g cm“^. The flight was terminated at 0710 CST on 
31 August and the payload recovered approximately 400 km west of 
Yorkton. 

2.2 Calibrations 

Both pre-flight and in-flight calibrations were carried out for the 
HEIST instrument. Most prominent in the pre -flight category was a run 
at the Lawrence Berkeley Laboratory heavy ion facility (Bevalac), from 
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which data were collected for the oonetruction of two~dimenelonal light 
collection efRciencj maps for DO through D7. Alio included in the 
pre-flight category were detailed electronic calibration! of the pulse 
height analyzers (PHA'i) used to convert PMT signals to digital form. 
In-flight calibrations have included cross -calibrations of the MWPC'i, 
residual light collection maps, PMT temperature-dependent nonlinearity 
corrections, and inter-detector normalizations. 

2.2.1 MWPC Calibrations 

The difference in arrival times of the delay line pulses, 6r, is 

6t 

linearly related to the position x: x * ~ 4*/3. The parameters a and /? 

vary from counter to counter and from sector to sector within a counter 
(pre-flight diagnostic tests indicated, for example, that a, whose 
nominal value is 14 nsec ^un~^ could vary by as much as 10% among the 
various counters). For a given MWPC our parametrization of response is 
6t— T_ 

r for the left sector 

«- 

x(<5t) s= for the center sector 

“o 

+r for the right sector, 

01 + 

where t_(+) = response of MWPC corresponding to the left (right) edge of 
the Csl stack, given a beam of normally incident particles; r = radius of 
Csl stack: To = response of MWPC to particle at center of MWPC; and 
ou,oc,a^.- delay of left, center, right sector, respectively. The 
quantities t_ and were determined for each MWPC from a Bevalac run 
(§2.2.3) of normally incident ^‘’Ar nuclei. The flrst-order estimate of the 
other parameters is a» =* a© * * 14nsec mm"^ etod tq = 0, for all 
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MWPC’i. 

We hare uaed thia firit-order eetimaie of MWPC response to 
determine, for each event, a particle trajectory a « m,z -f b,, 

y « m,z + b« by minimizing the quantities Xz *■ E and 

4 

Xr * S (yi ” ' where Zz(y)i* z-coordinate of MWPCX(Y)i. We 

iMl 

show in Figure 2.3a the distribution of the minimized Xy the flight 
data. The long tail for Xy> 7 mm is due to a contamination of the 
trajectory data by knock-on electrons, produced most likely In the 
0.69 gcm~^ aluminum shell which encloses the experiment, which 
confuse a delay line with more than one pulse. This background is 
discussed in detail in Appendix A. The peak at x^^ 2.7 mm indicates that 
the resolution for an individual uncontaminated y-measurement is 
CTy« 2.7mm (Bevington 1969). By adjusting the as and 
individually for each sector of each MWPC within the constraints imposed 
by counter fabrication (typical adjustments were a few percent, all were 
< 5/S) we have improved this distribution to that shown in Figure 2.3b, so 
that 7 jW imm. Shown in Figure 2.3c is the distribution of the 
minimized Xz» after adjustment of the parameters, for the x-trajectory 
measurement (the shape of the Xx distribution differs from that of Xj 
because we use only three measurements XI, 3,4 in determining the 
x-trajectory). The quantities Xz and Xy will be used later to reject events 
for which our determination of trajectory is inadequate for mass 
resolution. 
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Fifure 2.3 


MWPC reeolution. 

a) Xj distribution, first-order parametrization of MWPC response. 

b) Xj distribution, optimized parametrization 

c) Xk distribution, optimized parametrization. 


Counts/0.2 mm 
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2.2.2 PHA Calibraiions 

Prior to flight, electronic calibrations were :maue for the pulse 
height analyzers wnich are designed to logarithmically compress the PMT 
output. Diagnostic tests indicated a non-negligible temperature 
coedicient to PHA response, so calibrations were made at several 
temperatures and each PHA was assigned its own thermistor for 
recording of temperature in flight. Figure 2.4 shows typical variation in 
PHA output yrith respect to temperature given a fixed input voltage. A 
number of runs at different input voltages and temperatures provides us 
with a two-dimensional table of PHA response. For each event we have 
used linear interpolation in temperature and raw pulse height to convert 
PHA output for D0-D7 to a form which is proportional to the log of the 
PHA input, independent of temperature. 

2.2.3 Spatol Variations in Scintillator Response 

HEIST was taken to the Bevalac in February 1977. A defocused 
monoenergetic 900MeVamu“* penetrating ^Ar beam illuminated the 
instrument, allowing the construction of two-dimensional light 
collection efliciency maps for DO through D7. Baaed on the trajectory 
information supplied by the MWPC hodoscope, events were assigned to 
5 mm X 5 mm bins in x and y position at each scintillator. For each bin 
the mean of the pulse height distribution was calculated, and compared 
with the average pulse height over the entire crystal. The maps were 
then smoothed by performing a moving average over 25 mm x 25 mm 
areas to decrease the statistical error in the mean to 0.1— 0.2%. We show 
in Figure 2.5 the percentage deviation from average signal size as a 
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Figure 2.4 

Output of the D6 PHA vs. temperature for a fixed input signal 
(0.251333 volts). A 1 channel change in 'Pulse Height" corresponds 
approximately to a 0.126% change in input voltage. The figure indicates 
a temperature-dependent response of ~ -0.15%/”C. 


5260 
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ftinction of position on the crystal face for D6. As expected, Ught 
collection is most efficient near the edge of the crystal to which the light 
pipe is attached. If r(x,y) is the response at position (x,y), then the 
fractional rms variation in response over the face of the crystal is 

^ A “1]^ ] ® where r©*- ^J’dxdy r(x,y) and A» crystal area. 

_! 

The mean fractional gradient is - 7 “ fdxdy - 7 -^— r[(4^)^ + (4^)^]^- Both 

of these quantities are tabulated for D0-D7 in Table 2.1. Typical 
gradients are 0.1-0.2% n%m~‘ for the thick scintillators and 
0.2 -0.3% mm~* for the thin ones. 

A portion of this spatial variation in detector response is due to 
thickness variations over the face of the crystal and not to light 
collection variations. This rms variation (Table 2.1), however, is small 
compared to light collection variation, although it could in principle be 
deduced independently. While it is true that thickness variations ought 
tc be handled differently fx’om variations in light collection efficiency 
(the former should not be treated as a correction on measured 
scintillation light, but rather as an adjustment in thickness t in, e.g,. 
Equation 2.4), the error we make in lumping the two together is small 
(see § 3.5). 

Because the Bevalac calibration took place 1.5 years prior to the 
Yorkton flight, it might be expected that the response map measured 
there would be inadequate for application to flight data. Unlike the 
Bevalac calibration, however, there is of course no subset of flight data 
which consists of a particular nuclear species at a fixed energy and 
incident angle. The first step in using flight data for calibration 
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Figure 2.5 

Bevalac light collection efficiency map for D6. Nxxmbera represent the 
percentage deviation from the nominal response of the detector to 
relativistic ^Ar nuclei. Units are 2 % (i.e. a 'U" means between 0 and 2 % 
deviation, a "1" between 2 and 4%, etc.). The light pipe is attached to the 
crystal in the -x, +y quadrant. 
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purposes isolates a subset of that data consisting of a single charge (by 
methods described in Chapter 3). The data set consists of either 
particles which stop in a particular detector Dm. or penetrate all 
detectors D0—D7. In the former case, signal (denotes the 

PHA-cnlibrated output of detector Dm, corrected for spatial variation in 
response as measured at the Bevalac, equal to the log of the light output 
L„J is a good measure of particle energy. In the latte;' case we have no 
measurement of the particle’s total energy, only its energy loss in 
D0,Dl,...O7. One can take Equation 2.1 with the scaling in Equation 2.3 
to define the energy loss AE for a given M and Z in a particular detector 
as a function of the incident energy E. Differentiation gives 


aAEx 


s 1 - 


S(— ) 

S(“) 


[recall that S(c) is the specific ionization of a 


proton of energy c]. Thus AE is a measure of particle energy only to the 
extent that S(E’/M) is different from S(E/M). For penetrating events 
detector D7 best satisfies this criterion since, as the thickest detector, it 
gives the largest AE =* E-E*. Hence for penetrating calibration data sets 
we take signal D7 to be the best available measurement of particle 
energy. 


Given some calibration data set, then, we can remove variations in 
signal Dq correlated with energy and angle by fitting a function 
fn(seci),Dni) which predicts, for a given seciJ and Dm, the "nominal 
response" of detector Dn. We have used the form 

fn(seciy,Dm) * E E (seciS!)^"‘D^‘ ; 

l=ij=i 

the coefficients aj] are deter min ed by minimizing the quantity 
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[DB~fa(8eoi9,Dn)]^ (the sum i* over all events of the selected data 

mnU 

set). Once we have determined the functions f„, we can scatter plot 
Dq -> fQ(sec d.DoJ versus quantities such as position and temperature to 
look for correlations. 

As a specific example we show in Figure 2.6 a scatter plot of the 
deviation of signal D} from nominal response. D| -f|(sec‘d.D 4 ), for 
a-particles which stop in detector D4. versus the x* coordinate at 
detector Dl. (The x’-y’ axes are rotated 45® from the x-y axes; they 
correspond to the light collection symmetry axes of the scintillators. 
See Figure 2.2b.) The data have been corrected for light collection 
variation as measured at the Bevalac, yet a considerable spatial 
variation remains. This residual light collection variation is due most 
likely to a degradation in the surface polish of the scintillators over the 
~ 1.5 years between the Bevalac calibration and fiight. 

To remove such an effect for a given detector Dn, we fit another 
function 

<5n(x',y') « s s b,j (x')‘-»(y')^-» 

i=ij=i 

whose coeflficients by are determined by minimizing the quantity 

E [Dn-fn(sec'«5'.Dm)-'<5n(x', y')]®. where f^ has already been determined, 
eventa 

The quantity is then the multiplicative correction factor which, 

when applied to Lb = c°" for an event at (x’.y’) on detector Dn, will remove 
the residual spatial dependence of output Lq. 

We have used this technique to look for residual spatial variations 
in the response of detectors D0-D6. For D4-D6 we have used 
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Fifure 2.8 

Deviation from nominal response of detector D1 vs. x* for a-particles 
stopping in detector D4. The D4 light pipe is in the -hx' direction. Since 
'Dj" is a logarithmic measure of the output of Dl, full scale on the 
ordinate corresponds to a factor of 1.8. 


I 


3 
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penetrating oxygen data for the calibration (n»4,S,6; m«7) under the 
assumption that the residual variation in D7 is negligible (the 
consequences of this assumption are discussed in §3.5). Fo.‘ D0-D3 we 
have used a-particles stopping in detector D4 (n«0.1,2,3; m«4), having 
corrected D 4 for residual variation with the oxygon data set. We list 
under '’residual" In Table 2.1 the fractional rms variation s„ in residual 

response for D 0 -D 6 , where s^ ■ fdx' dy* [ - — 1 ]® with 

A r Q 

r'(x’,y’)* and r'o ■ ^/dx'dy'r'(x',y'). Also tabulated are the 

X 

mean residual gradients fdx' dy' J f fV • Finally, 

these same quantities corresponding to the superposition of the Bcvalac 
and residual response maps are tabulated under the column in Table 2.1 
labeled 'Total". Table 2.1 thus summarizes the spatial nonuniformity 
characteristics of our scintillators. 

The precision with which we can make the residual correction 
depends on the statistical accuracy in the determination of the 
coefficients by for a particular (Sq. That is, the quantity J)[ Dn-fn-dn 
is not precisely zero, because other effects (for example, Landau 
fluctuations; see §3.5) cause scatter in signal Dq. Although the exact 

value varies with (x',y’), typically the correction is known to ± Sn (— j 

where N is the number of events used in determining dj,. This value 
ranges from 0.2-0. 3% for D 0 -D 6 . 


1 


S3 

2.2.4 Tamperature-Depmadtrit PMT Output 

There ars two effecti which can cauie PMT output to change with 
respect to temperature given a fixed amount of energy deposited in the 
aesociated scintillator. One is the temperature dependence of the 
scintillation efficiency of Csl which is approximately - 0 . 4 X/‘*C (Birks 
1964). The other is the temperature-dependent gain of a PMT which is 
expected to be of the same order (Wiedenbeck 1978). PMT temperatures 
were monitored in flight and ranged from 30 to 40 *C, The Csl 
temperature was not measured directly, but nearby measurements 
indicate that its total variation in temperature was similar to that of the 
PMT's. Based on carbon stopping in D7 we have derived temperature 
, coefficients for DO through D 6 (Table 2.2), typically ~-l%/®C, under the 

I 

I assumption that the temperature coefficient for D7 is the average of 

those for D 0 -D 6 (an iterative approach is thus required). The technique 

i 

! was identical to that used in determining the residual spatial variations 

in scintillator response, with temperature replacing position as the 
^ independent variable. The variances of the fits and number of events 

I used in the calibration have limited the precision with which we can 

correct PMT output for temperature variations to ± 0.4%. 

2.2.5 Inter-Detector Normalization 

' Because we add signals from different detectors together for 

analysis of mass (e.g. as in Equation 2.4 where we add the separate 
measurements E’ and AE to get E), it is necessary that the detectors be 

( normalized among themselves. We show in Figure 2.7 a histogram of 

signal D 5 , normalized to vertical incidence, for penetrating particles. 


IL 


33 


Table Z.2 PBfT Tcmperatiure Coefficients anit Detector Normalizations 


Detector 

Temperature Coefficient 

Normalization 

DO 

-0.56 V“C 

1.53 

D1 

-1.23 

0.91 

D2 

-0.76 

0.75 

D3 

-0.65 

1.03 

D4 

-1.05 

0.93 

D5 

-I.IB 

1.16 

D6 

-0.94 

0.37 

D7 

-0.91 

1.20 


\ 
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Ficure 2.7 

Penetrating distribution of signal D 5 = InLs (the units of L are defined in 
the text). The cutoff at D 5 W 2.7 corresponds to the selection for this plot 
that the high level discriminator not fire. 
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[Dq b inLai the units of L, which we call "normalized light units" or nlu, 
are such that 1 nlu corresponds to a charge signal at PMT dynode 7 of 
nominal magnitude llpC (IpC* 10"^^ Coulomb)]. The peaks correspond 
to boron, carbon, and oxygen (nitrogen, about one-quarter as abundant 
as carbon, is not resolved). The major contribution to the width of the 
peaks is the variation in particle energy at the top of DO. The most likely 
pulse height for a given charge, however, corresponds to minimum 
ionizing /? (the minimum dE/dx in Equation 2.2 which occurs at 0.95 
or E/M w 2 GeV amu"^) and is therefore proportional to detector 
thickness. The unnormalized PMT output of detector Dn due to a 
minimum ionizing particle of a given charge which creates an amount of 
light per unit distance 1 is qn = gn^^^« '»^here g„ = gain of PMTn, 
Cn - fraction of light created which is transported to the PMT face and 
converted to photoelectrons, and t^ = thickness of detector Dn. The 

Cjj q„ 

normalization constants c„ are chosen so that — - is the same for all 

8n ^ 

detectors: they are tabulated in Table 2.2. For those detectors for which 
we have three independent determinations of c„ corresponding to the 
boron, carbon, and oxygen peaks, the Cq's agree to better than ± 0.01. 
Given our method of determining the absolute light to energy 
relationship (§3.3), this error contributes to errors in mass only in 
second order, and is therefore negligibly small. 
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Chapter 3 

Maaa Analyais and Resolution 

Up to this point we have been concerned primarily with the various 
corrections which must be made to the raw data to put it in a form 
suitable for analyzing masses. In this chapter we deal with a more 
fundamental problem, namely the relationship between the amount of 
energy deposited in Csl and the corresponding scintillation light, and the 
dependence of that relationship on particle properties such as charge, 
mass, velocity, and specific ionization. Having solved this problem we will 
be able to deduce the isotopic compositions of boron, carbon, and 
nitrogen as measured by our instrument at 5gcm”^ atmospheric depth. 
Finally we discuss those factors which contribute to the mass resolution, 
and compare the expected mass resolution with the performance of the 
instrument. 

3.1 Theory 

We begin by considering in more detail the AE—E' technique on 
which our experiment is based. If R(E,M,Z) is the range in Csl of a 
particle of total kinetic energy E, mass M. and nuclear charge Z, then 
after traversing a thickness t at angle tJ in which the particle loses an 
energy AE, the residual range is R(E’,M,Z), where E' = E— AE (Figure 2.1 J. 
The range-energy function of a heavy ion is related to that of a proton 

(Equation 2.3) by R(E,M,Z) * R_(^). We thus have 

D / AE + E\ _ /E\ _ Z^tseciJ /o i\ 




We show in Figure 3.1 this theoretical relationship between AE and 
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E’ for various particles which stop in D7 of our instrument. We have 
taken t»13mm (the thickness of D6), and used the tabulation by Janni 
(1966) of the range-energy function for protons in Csl. AE in 
Equation 3.1 corresponds to the energy lost in D6, and E’ corresponds to 
the remaining energy lost in 07. We have assumed (solid lines) and 
have plotted one 'track" with i>s30° (dashed line) to show the angular 
dependence. Because of the T? scaling of dE/dx (Equation 3.2) we see 
from this figure that the separation between tracks of adjacent elements 
is much greater than that between isotopes of a single element. One can 
calculate from Equation 2.4, which is based on the power law 
approx.imation Rp(c) =kc*, that 

9AE N _ a-1 Z 6AE \ 

aM " a+1 M aZ E',^’ 


With a>« 1.66 for Csl (Wiedenbeck 1970) and M = 2Z this yields 

•^^)e*z = 0.12-^^) 7 . 

aM az 


An exact calculation shows that the separation in AE between and 
at E' = 1 GeV is 6.4% of the separation between and at the same E'. 
Thus a change in one mass unit is equivalent to a change of 0.128 charge 
units. 


3.2 Charge Identification 

We now describe the means by which we can identify in our data 
the charge of each particle which stops in D7. Figure 3.2a is a scatter 
plot of Lc vs. hf for particles stopping in D7. We select particles that stop 
in D7 by requiring that neither D8 nor the anticoincidence detector be 
triggered, as well as requiring that the D7 pulse height trigger the high 
level discriminator. No corrections have been made to the data beyond 
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Figure 3.1 

AE = Ee vs. £'*£7 for boron, carbon, nitrogen, and oxygen isotopes 
stopping in detector D7 with an incident angle 0 ° (solid lines). Also 
shown is the 'track" for at an incident angle of 30“ (dashed line). 


2000 



41 


Fiffure 3.2 

Le vs. Lr for particles stopping in detector D7. 

a) Raw data (except for PHA calibrations, and inter-detector 
normalizations). 

b) Light collection maps, PMT gain shifts, seciJ correction. 


/sec® (niu) 1-6 (niu) 
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the inter-detector normalizations aKid PHA calibrations discussed in 
Chapter Z. Even at this stage one can see poorly resolved boron, carbon, 
and oxygen element tracks. The data In Figure 3.2b have been corrected 
according to the Bevalac and residual light collection maps on the basis 
of their trajectories, and corrected for PMT gain changes on the basis of 
the appropriate PMT temperatures. In addition, a first order pathlength 
correction has been made to by normalizing it to vertical incidence, 
that is, dividing L« by seed. This correction ignores higher order effects 
due to the fact that the particle is slowing down as it traverses D6 so 
that its specific ionization is increasing. It is, however, sufficient for 
charge resolution. One can see in Figure 3.2b clear element tracks 
corresponding to boron, carbon, nit»-ogen, and oxygen. 

We have assigned a charge to each event in Figure 3.2b by 
empirically determining a family of curves which parametrize the 
tracks: 

InZ = Aln-^„+B InL' + C , 
seed 

where A, B, and C are independent of Z. The form is suggested by 
Equation 2.4 with M = 2Z and AE<<E'. In addition to assigning a charge 

to each event based on its location in the - - -L7 plane, we can plot 

sec V 

vs. L« + L7, vs. L4 + L8 + L7, and so on (we have not included L5 

seev secij 

because for nitrogen and oxygen the signal out of D5 is between gain 
ranges of its PHA). For each of these plots we have also empirically 
determined a family of curves which parametrize the tracks, so that we 
have six independent charge assignments Zj, 5?* i=0 , for each event. 
We are thus able to eliminate events which undergo a charge-changing 
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Xt 

Quclsar interaction in the CsL In Figure 3.3 we show the quantity 

"art 

plotted against where *uad xl" 2^i(Zi“Z»y,)® (the 

weights wi are determined from the variances of the fits to the tracks, 
and have been normalized so that ][)w(* 1). There is a clear depletion of 
non~integral events as well as a general improvement in charge 

resolution for small ■= — . That the means of the four charge groups do 

^■rt 

not all occur at exactly integral values of is due to the 

approximations inherent in the choice of the form used for the 
functional dependence of Z on AL and L’ (in particular the 
approximation that A. B. and C are independent of Z). 

We show in Figure 3.4a the subset of events in Figure 3.2b for 

Xz 

which 4.5< Zav«< 0.5 and •=r^< 3.5%, together with the nominal tracks for 

^ave 

Z = 5, 6, 7,0. Figures 3.4b-f consist of the same events as in Figure 3.4a, 

and show — vs. 2 i=0,l,...4. Charge histograms 

seen? 

corresponding to Figures 3.2a, 3.2b, and 3.4a are shown in 
Figures 3.5a-c. Based on the total inelastic cross section for particles in 
Csl we expect that 20-25% of the events with 5^ Zi 8 will undergo a 
charge-changing nuclear interaction in the Csl before stopping in D7, 
This is consistent with the attenuation in events from Figure 3.5b to 
Figure 3.5c. 

From Figure 3.5c we also see that the typical charge resolution is 
ozwO.17 charge units. From our earlier discussion this is equivalent to 
a mass resolution of only 1,3 mass units. Although we have 
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Fifure 3.3 


— — vs. The boxed regions define a selection criterion for mass 

^avB 

analysis (§3.4, Table 3.1). 


6 



j 


ove 



47 


Fifure 3.4 


— vs. 2 Li for 4.5<Z,y*<8.5 and 3.5!J. Also plotted are the 

seClJ 5#J>1 ^ave 

nominal integral charge tracks. All L’s are in nlu, 

a) is0. 

b) is4. 

c) is3. 

d) i»2. 

e) i»l. 

f) i*0. 
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Figure 3.5 

Charge histograms (Zs) for particles stopping in detector 07. 

a) PHA calibrations and inter-detector normalizations only; corresponds 
to Figure 3.2a. 

b) Corrections for spatial variations in detector response and PMT 
temperature, in addition to above. A first-order correction for 
pathlength has been made by dividing Lq by sectf; corresponds to 
Figure 3.2b. 

c) Same ns b), except require that 4.5 < Z*v«< 0-5, and ~ — < 3.5?;, The 

attenuation in the 'peaks'* (defined to be the four charge bins for a given 
element with the most counts), from b) to c) is 26± 2%; for the 'Valleys" it 
is 45± 3%; the total attenuation is 31 x 2%. 
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demonstrated that we can identify the charge of each particle on an 
event -by-event basis, our analysis must be more sophisticated if we 
hope to resolve isotopes of a single element. We note here that we have 
not at this point made any selection of events on the basis of the quality 
of trajectory data because our sec correction is sufficiently naive that 
to do so would not improve the charge resolution significantly. 

3.3 Scintillator Saturation 

The ideal situation would be one in which our detectors were 
perfectly linear, that is, the scintOlation light would be directly 
proportional to the energy deposited. Then a single constant would 
suffice to correct our measured signal to an energy, and we could simply 
use Equation 3.1 to calculate the mass for each event. The well known 
phenomenon of scintillator saturation (Birks 1964) means, of course, 
that this is not the case. From Figure 3.1 one can calculate that, at 
E' = 1 GeV, AE for oxygen is a factor of 2.8 greater than that for boron. 
On the other hand, in. Figure 3.4a at L7 = 30 nlu, for oxygen is only a 
factor of 2.2 above that for boron. This is already a first order 
measurement of heavy ion saturation in Csl. 

To deduce the light-energy relationship for Csl we do not assume 
that such a relationship is necessarily the same for different elements 
because a particle with a given charge and total kinetic energy has a 
"Bragg curve" (dE/dx vs. penetration depth) very different from that of a 
particle with the same energy but a different charge, and we expect that 
the differential scintillation efficiency dL/dE will depend mainly on 
dE/dx. We thus coniine ourselves to a single element at a time, and use 
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the carbon data to illustrate our procedure. Since we expect that 
carbon will consist primarily of we assign to each event with Z»6 a 

mass Ma*12amu. Our goal is some functional relationship -p s a(L), 

where L (in ’hlu"; see §2.2.5) is the total light out of a scintillator in which 
a particle of total kinetic energy E (in MeV) deposits all of Its energy. 
Our zeroth order (false) assumption is that a(L) » constant » oq (perfect 
linearity). Then Equation 3.1 becomes 

00 ( 4 ^ ] - (3 8 )- 

Having set M~12amu, all quantities in Equation 3.2 are known except 
oq. For each event, then, we calculate the value of a© required to satisfy 
the equation. We then scatter plot Oo vs. L’ and AL + L' (Figure 3.6a) and 
see that our assumption of constant a is false (as expected). However, 
we can use the scatter plot to define a function (5oi(Li,L2) a ^(Lj) - o’(L 2 ) 
by fitting a straight line oT to the oq vs. L' portion of the plot. The 
function da teUs us how different E/L is for Lj = L' + AL and La = L'. 


For our next iteration we no longer assume, then, that a is 


constant, and Equation 3.2 becomes 




I -t- L') 1 ^ D r g(L*) L* -t _ Z^tsecT? 


(3.3). 


Again, every quantity in Equation 3.3 is known for each event (we derived 
da from our first iteration) except a, so we again calculate a event by 
event and scatter plot this vs. L'. We continue iterating in this fashion 
until the da function we put in is consistent with the cT function we get 
out, in which case cr(L) is E/L for In Figure 3.6b we show a scatter 
plot of a(L') and a(AL + L') vs. L’ and AL + L' as deduced from the last 
iteration. In addition we show the straight line fit or(L) deduced from the 
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Figure 3.6 


Determination of a(L) = — for carbon. Crosses correspond to L'^L?, 

circles to AL+L' = Ls+L^+Ly. 

a) First iteration (assumes a = constant), 

b) Convergence. 

The selection criteria are 5.7 6.3, secd< 1.2, and Xx.y< 5 mm. 
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previous iteration from which came 6a used in the last iteration. The 
agreement is good, indicating that the procedure has converged. 

For boron and nitrogen the procedure for determining the E/L 
fimction is identical, except that, since these elements are expected to 
consist of substantial quantities of more than one isoto pe, we assign to 
each event of the given element the same average mass arbitrarily. 
Since we expect boron to consist of and we assign 10.5 amu to 
each event with Z^S. Simiiarly, nitrogen consists of and so we 
assign 14.5 amu to each event with Zs . ‘ vfiU take into account this 
arbitrary assignment of average mas. in determing the E/L 

function by letting the offset of the mass scale be a free parameter later 
on. We show in Figures 3.7a, b the convergence of the iterations for these 
elements. 

We have found that there is a small but non -negligible dependence 
of a on angle d, and have taken this into consideration by adding another 
term to the form used for Specifically, we have used the form 
^■(L.d) = AL + B+Csecd, and find that C> 0. As seed varies from 1 to 1.2, 
o’ varies by typically 3%. Several candidates have been identified which 
might explain this dependence, including i) systematic overestimate of 
the nsec mm“* delay in the MWPC delay lines (§2.2.1), ii) thin dead layer 
in the AE-detector, iii) errors in inter-detector normalizations (§2.2.5), 

and iv) systematic error due to the assumption that — is related linearly 

Lj 

to L (see below), none of which is in itself sufficient to explain the effect, 
although all have the right sign. 

In Figure 3.0 we show the relationship between dL/dE as deduced 


N. 
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Fifure 3.7 

E 

Determination of a(L) » y for a) boron and b) nitrogen. For both 

Lf 

elements, L' (crosses) = hj. Circles correspond to AL + L'. For boron, 
ALssLg + Lo, for nitrogen, AL = L«. For boron we require 4.6 < Z,y,< 5.4, for 
nitrogen 6.8 < 7.4. For both elements only events with seciS< 1,2 

and Xx,j < 5 mm are used. 
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from our E/L functions and the specific ionization dE/dx for boron, 

carbon, and nitrogen, using 

dL.,dEj»,. 1 . 1 

dE 'dL' d/. Es 2AL + B+Csecil ' 
dL^^L^ 

The curves have been normalized to dL/dE for minimum ionizing carbon 
(Figure 2,7). Given that the size of the error bars (which depend on the 
statistical errors in A, B, and C) are greater than or of the order of the 
separation between charges, our data are consistent with dL/d£ 
depending only on dE/dx, independent of charge and detailed structure 
of the distribution in energy transfer to electrons in the ionization 
energy loss process. 

The shapes of the curves in Figure 3.0 are a consequence of the 
E 

assumption that ~ is related linearly to L, and should not be interpreted 

too literally [other investigations (see Birks 1964 for a summary) suggest 

that ~ might be linear in ln'~]. To estimate the magnitude of the 

error we make in mass due to this assumption we have performed an 

analytical calculation which assumes i) ® a + b In 4 ^ (specifically, a 

dE dx 

straight line on Figure 3.0 which connects the endpoints of the carbon 
curve) and ii) the power law approximation Rp(c)=ke*. With these 
assumptions, we calculate AL and L’ pairs, ror a given nuclide, as a 
ftmction of the range R 7 in D7. These pairs are then analyzed, just as 

was the flight data, to derive a straight-line ~ vs, L relationship 

Jut 

P 

[because of assumption i), -jj vs. L is in fact nonlinear]. We then use the 
derived linear relationship to convert the L's into E’s, and finally 


1 . 
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Figure 3.B 

Differential scintillation efficiency dL/ dE vs, specific ionization dE/^*x for 
boron, carbon, and nitrogen. The dL/dE axis has been normalized to 
minimum ionizing carbon (Figure 2,7), where dE/dx*20,3 MeV mm"**. 
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calculate masses using Equation 2.4. We find that the resulting rms 
variation in mass, for 2mm^ 17mm and secn9$ 1.2, is O.OSamu 
for 0.07 amu for and 0.1 1 amu for 

Given an E/L function for an element, hovr might we expect it to 
differ for different isotopes of that element? A reasonable assumption is 
that dL/dE depends only on dE/dx, since two particles with the same 
charge and specific ionization (and therefore the same velocity) give rise 
to identical electron energy transfer spectra. Since the velocity is 
equivalent to kinetic energy per unit mass, we have 

§(E.M,Z) = tz(|). 

SO that 

E . E/M 

L(E,M,Z) = / dE-^=M/ dctzW 

no aL •'o 

upon change of variable c = E/M. For another isotope of mass M' we have 

E/M* 

L(E,M\Z) = M' r de fz(c), so that 

0 

L(E.W.Z) =^L(E^,i?,M.Z) (3.4) 

where we now explicitly include the iJ-dependence discussed earlier. 
Thus once we have our E/L function for the average mass of an element 
we can use Equation 3.4 to tell us what it is for particular isotopes. 

3.4 Mass Analysis 

We can now calculate masses on an event by i basis. The 
procedure is to calculate M such 

(3.5). 

[we use E(L,tJ,M,Z) here to denote the total kinetic energy E o.f a particle 
of mass M and charge Z at incident angle i? which creates a total amount 
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of light L as it comes to rest in a scintillator], using our E/L function and 
the scaling implied by Equation 3.4 to calc .: . E(L,i9,M,Z) given a 

measured L and -A The problem with this approach is that it uses 
unphysical nonintegral masses in the conversion from light to energy. 
Of course we have no a priori knowledge of the mass to tell us which is 
the correct M to use in the conversion. Since our ultimate goal is to 
obtain relative abundances of the various is' opes, we can incorporate 
this lack of knowledge directly into our maximum likelihood analysis of 
the mass distribution, which we now describe. 

Suppose a given element consists of two isotopes, fraction fi of 
mass mj and fraction 1-fi of mass m 2 . We assume our resolution 
function for a single isotope is a gaussian with variance For each 
event we calculate, according to Equation 3.5, two masses, M} assuming 
mj in the L to E conversion and M 2 assuming m 2 in the L to E conversion. 
The probability that a mass distribution characterized by a^andfi 
results in the measurement Mj and M 2 is then 

p(Mj,M 2 ) = ^ ' [fiexp[- — ^ ] ■+• (l-fi)exp[- ^ ^ ^ ■ ] ], 

vZn a 2a^ 

A maximum likelihood analysis consists of forming the likelihood 

function L(fj) = f] pj, where pj is the probability of the jth event. The 

maximum likelihood estimate of f| is that which maximizes L(fj). The 
extension to more than two isotopes is straightforward. 

We have used this technique to calculate maximum likelihood 
estimates of the abundances ^®C/C, and ^®N/N. Because the gain 

of our mass scale is fixed by the mass scaling of the light to energy 
function, we have not allowed mi-m 2 to vary in the analysis, but to take 


\ 


65 


into account our arbitrary assignment of average mass in determining 
the E/L functions, we take the offset of oxir mass scale to be a free 
parameter. In addition we allow o to vary, so that the likelihood 
function depends on fi.mi.andcr, with fixed mj-ma. For carbon we 
expect roughly 1% of the events to be due to nuclear interactions in 

the atmosphere above our instrument (§4.2), so that for this element we 

Ur 

have allowed the quantity --7- to vary between 0 and 2%. 

For boron and -carbon we have used the sum of L5 and as the AE 
measurement and Ly as the E’ measurement. For these elements we 
have also separately calculated masses using i) D6 as the AE detector 
and D7 as the E' detector and 11) D5 as the AE detector and the sum of 
D6 and D7 as the E' detector, allowing us to reject events if these masses 
are inconsistent. Because we do not have L5 for nitrogen, we have only 
the mass i) in that case. 

We have made scatter plots of MU3 (the notation indicates the 
mass calculated using Equation 3.5 assuming the subscripted isotope in 
the L to E conversion), Miag, and Misjj versus such quantities as seciJ, L5, 
Lfl, and L7 to be sure that our mass distributions are not biased with 
respect to these quantities. In addition, we have plotted mass versus 
Xx and Xy to see whether our resolution can be improved by rejecting 
events with large values of these parameters which characterize the 
quality of trajectory data (§2.2.1, Appendix A). 

We show in Figures <3. 9a-d scatter plots of seciJ, AM [denoting the 
difference between the masses i) and ii) above], Xx* Xy versus Mu^* 
and similar plots for M113 and M»jj in Figures 3.l0a-d and 3.11a-d, 
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respectively. Because we have reasonable statistics for carbon, our 

selection criteria for all elements are to a large degree based on the 

carbon data. The cuts at seci9sl.2, AM = ±0.Bamu, and 

occur at or near values where there is a discontinuity in the distribution 

of the quantity. There are two cuts for Xy» depending on whether three 

or four y-MWPC measurements were xised in determining the 

y-trajectory (no events appear with Xy> 3 mm when four 

y-measurements were used because such events were re -analyzed after 

discarding one of the four; see Appendix A). For boron we have excluded 

events for which Xy> 2 mm when only three y-measurements were used 

because two of the three events so excluded (Figure 3. lOd) have 

Miig< lOamu, and the exclusion significantly improves the error in the 

10b xz 

determination. The slightly smaller value of the ■= — cut for boron 
o ^ave 

eliminates a single event with an anomalously low mass. The Z„e outs 

are all of width 0,6 cheurge units and are centered approximately on the 

mean Z*ve t-he particular element (Figure 3.3). From Figure 3.1 lb we 

see that most of the nitrogen events with an anomalously low or high 

E 

mass are correlated with small hf, a region where the — function is 

poorly known (Figure 3.7b), We have thus excluded nitrogen events for 
which Li 7 < 20. 

In Table 3,1 we summarize the selection criteria, parameters used 
and range of variation, and results of our maximum likelihood analyses, 
and in Figure 3.12 show mass histograms based on those results. Each 
histogram has two mass scales associated with it corresponding to which 
isotope is assumed in the conversion from light to energy. (We can do 
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Fieur« 3.9 

Carbon selection. Scatter plots of a) seci?, b) AM, c) d) Xy vs. Mit^. 

Circles in d) indicate events for irhich all four y’s were used in 
determining the y-trajectory; crosses indicate that only three y's were 
used. The selection criteria are indicated by the hatched lines. 




69 


Figure 3.10 

Boron selection. Scatter plots of a) seciJ, b) AM, c) x*. and d) Xy vs. Mug. 
Circles in d) indicate events for which all four y’s were used in 
determining the y-trajectory; crosses indicate that only three y’s were 
used. The selection criteria are indicated by the hatched lines. In d), 
the hatched line applies only to the crosses. 
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Fifurte 3.11 

Nitrogen selection. Scatter plots of a) seciJ, b)Ljr, c) Xz> d) Xy. vs* Misj^, 
Circles in d) indicate events for which all four y’s were used in 
determining the y-trajectory; crosses indicate that only three y’s were 
used. The selection criteria are indicated by the hatched lines. 
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Figure 3.12 

Mass histograms for a)carbon, b)b?ron, aad c)nitrogen. The different 
scales for the horizontal axes are labeled according to the assumed mass 
in the L to E conversion, as explained in the text. 
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this because the difference between Mtsg and Mug, for example, is 
constant to within ± ~ 0,1 amu.) The statistics and resolution for 
nitrogen are insufficient for the maximum likelihood method to deduce 
the offset of the mass scale, so we have not allowed m^ to vary for that 
element. We thus rely on the fact that both theoretical calculations and 
measurements by other investigators (§4.3) indicate that our assignment 
of 14.5 amu as the mean mass in the determination of the nitrogen E/L 
function is correct. The masses for carbon in Figure 3.12 have been 
adjusted by only 0.02 pmu according to the maximum likelihood 
estimate of mj; for boron there is no adjustment because the maximum 
likelihood estimate of mj is equal to the actual mass of namely 
9.94 pmu. 

From Figure 3.12 we see that, although the statistics are limited, 
the mass resolution at boron is sufficient to validate our assumption of 
how the E/L function scales for different isotopes of the same element, 
because the two peaks occur at the right mass values. That is, the lower 
mass peak occurs at 10 amu on the mass scale and the higher mass 
peak occurs at 1 1 amu on the mass scale. We also see in this figure 
that o'lr mass resolution, and consequently the precision with which we 
can determine isotopic abundances, is degraded considerably as we go 
from boron to nitrogen. 

3.5 Mass Resolution 

The factors which contribute to the mass resolution in our 
instrument can be divided into two categories. The first consists of 
those processes which place a fundamenLal limit on the AE-E' technique 
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of resolving isotopes in that they cannot be reduced by increasing the 
precision of our measurements. Multiple Coulomb scattering (mcs) and 
statistical fluctuations in the ionization energy loss process ("Landau 
fluctuations") comprise this category. The second category consists of 
uncertainties in our measurements of particle trajectory and 
scintillation light. Decreasing these measurement errors will increase 
our mass resolution. 

To evaluate the contribution of the various factors to the mass 
resolution we have used a Monte-Carlo technique which simulates 
multiple Coulomb scattering, Landau fluctuations, measurement of 
particle trajectory, spatial variations in light collection, and other 
miscellaneous effects. Events are generated with whichever of these 
processes we wish to include, and then analyzed for mass according to 
Equation 3.1. From the widths of the resulting mass distributions we 
determine the particular contribution to the mass resolution. We have 
taken D5+D6 as the dE-detector for boron and carbon, and D6 as the 
AE-detector for nitrogen (as in the analysis of flight data). For all 
elements we have taken D7 as the E’-detector. The event selection in 
the simulation is essentially identical to that used in the analysis of the 
flight data. The results are given in Table 3.2. A discussion of individual 
processes together with illustrative calculations baicd on Equation 2.4 
follows. 

In traversing a finite thickness ox matter a charged particle 
undergoes many small deflections as .it interacts with electrons in the 
matter. This multiple Coulomb scattering (Rossi 1952, Bischei 1972, 
Jackson 1975) means that a particle impinging on a detector at incident 
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Table 3.2 Contributiona to IfaaB Resolution 



tOQ 

‘*C 

‘^N 

Contributor 

Multiple Coulomb 

0.13 amu 

0.14 

0.10 

Scattering 

Landau Fluctuations 

O.OB 

0.09 

0.12 

Measured seed 

0.02 

0.02 

0.02 


(0.23)t 

(0.20) 

(0.34) 

Light Collection 

0.0b 

0.06 

0.00 

(0.10) 

(0.21) 

(0.27) 

Photoelectron Statistics 

0.05 

0.06 

0.07 

Calibration Statistics 

0.09 

0.09 

0.14 

Monte Carlo Total 

0.20 

0.21 

0.20 


(0.32) 

(0.30) 

(0.46) 

E/L linear in L 

0.05 

0.07 

0.11 

No D7 Residual Map 

0.16 

0.19 

0.20 

Total 

0.20 

0.29 

0.41 


(0.36) 

(0.42) 

(0.55) 

Achieved^ 

0.32 

0.37 

0.50 


± 0.06 

± 0.04 

± 0.10 


fValues in parentheses were calculated assuming a degraded resolution in 
MWPC XI of tr = 30 mm (Appendix A). 


jFrom Table 3.1. 
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angle ^ will generally emerge at some slightly different angle i9*. The 
mean square angle of scattering is given by (Bischel 1972) 

<t| = < (U'-0»> = 1.7BX io-7r»d»cm»g->-?%2al^£l (3.6) 

where 


Z = atomic number of parade, 

M = mass of particle (pmu), 
s particle velocity (in units of the speed of light), 

7 = (1 

Zm- Aro = mean atomic number, freight of detector material, 
p = density of material (4.51 g cm"^for Csl), 

X » thickness of material (cin), 

B = weak function of Z^, Z, pz. 

For particles and energies of interest in our experiment B is Independent 
of /3 and is given by Bc«i = 2.5 logio(px cm^g”^) + 12.4. Equation 3.6 
applies if the thickness of the material is such that the particle velocity 
does not change significantly in traversing the material, generally not 
the case in our experiment. We have derived (Appendix B) a formula 
based on Equation 3.6 which includes the fact that /? is not constant; for 
lOgcm”^ this represents about a 50% increase over Equation 3.6, 
because the term 7 ^^ in the denominator decreases as the particle 
slows down in the material. Scattering by an amount o^in the Csl above 
the AE-detector results in a pathlength uncertainty in that detector of 

seci3 ■■“cr”’ 


whereas scattering by an amount in the 


AE-detector results in a pathlength uncertainty of 


g,eei» _ tan 
seci3 ”~v¥* 


From Equation 2.4 we have 

_ 1 9M \ _ _ 1 

'■“M - M 3seci>^®*'^^ a-1 seciS* 

T 3 rpical rms scattering angles re cr^w25mrad in DO— D4, and 
cr^RsSOmrad in D5-D6, giving 0. 14 amu for M=12amu and t?=20°, 


BO 


consistent with the Monte-Carlo calculation. [Because D6 (and not 
D5+D8) was used as the AE-detector for nitrogen, entries in Table 3.1 
for nitrogen are generally higher than those for boron or carbon.} 


The second fundamental limit to our mass resolution is the 
statistical nature of the ionuation energy loss process. These Landau 
fluctuations mean that a beam of monoenergetic particles of a given 
charge and mass incident on a thin detector will emerge with a 
distribution of energies, that is. they do not all lose the same amount of 
energy in the material. For particles and energies of interest in our 
experiment the distribution in energy deposited AE is a gaussian with 
variance (Jackson 1975) 

- (0.396 MeV)2cm2g-‘ ^ pxy^l - -^)D^ (3.7) 


where all quantities are defined as in Equation 3.6 except the 
"deceleration factor" D which takes into account the slowing down of the 
particle if x is sufllciently thick so that AE Is not small compared to the 
particle's total energy. We have used a form for D derived by Wiedenbeck 
(1973) which ranges from unity at x/R = 0 (R is the total range of the 
particle) to 2 at x/Rs=0.9. Landau fluctuations have the effect of moving 
events which are on a nominal isotope track (such as in Figure 3.1) off 
that track along a line of constant AE+E'. Equation 2.4 gives 


3M 


— \ „ _ a iw 

Landau “ /E.aeed ~ g 


(l-iL) a 
M ^ R'' 


_x 

R 




A 200 MeV amu~* nucleus incident on D5 has a range in Csl of 
34.5 mm. With D5+D6 as the AE— detector x/Rw 0.58 for which D«1.2. 
From Equation 3.7 we calculate cr^E = 6.7MeV so that au.ijjidau'® 0-10 amu, 
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again in agreement with the Monte*>Carlo calculation. 


Errors associated with the finite precision with which we measure 
particle trajectory cause an imcertainty in the particle's 1) angle of 
incidence, which in turn results in an additional (beyond multiple 
Coulomb scattering) uncertainty in the pathlength in the AE detector, 
and ii) absolute position at a scintillator, resulting in an uncerteunty in 
the correction of scintillator output due to spatial variations in light 
collection. Contribution i) depends only on the spatial resolution of 
individual MWPC's, and for our geometry is given by 


= 0.0016 mm“* sinHiJ ounc (Wiedenbeck 1978). For 

sec V 

T>= 20®, and M = 12 amu we have 

‘^kiect»= “ 0.02 amu. 

Contribution ii) depends additionally on the magnitude of the spatial 
gradients of the light collection variations (Table 2.1): 

3M V 

O’AE 


1 


= M 


l-(l-^)‘ „ 
a ' R f^AE 


a— 1 


2 ^ 

R 


AE 


(3.8). 


With D54-D6 as the AE-detector, 0.15% mm"^ gradients, 1mm MWPC 
resolution, 12 amu, and z/R=0.6 we have 0.03 amu. 


In our Monte-Carlo analysis we have evaluated these contributions 
to the mass resolution due to errors in trajectory measurement with two 
different assumptions of MWPC resolution. One, we assume all MWPC’s 
measure particle position accurate to the nominal ± 1 mm (the result is 
consistent with the rough calculation above). Two, we assume the 
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resolution of XI is degraded, because of knock’-on electrons, to ± 30 nun. 
with all other MWPC's at the nominal value. The choice of 30 mm is 
somewhat arbitrary, as we do not know precisely what the effect of 
knock-on electrons will be on MWPC resolution for various particles and 
energies, but 30 mm is suggested based on our analyses of the 
y-trajectories (Appendix A). We choose XI as the 'bad" coordinate since 
events are not analyzed if X3 and X4 are inconsistent (X3 and X4 are 
separated in the vertical dimension by only 4.5 cm whereas we have no 
x-measurement at the top of our instrument other than XI). In 
Table 3.2 we have put in parentheses those contributions to mass 
resolution which have been evaluated assuming 'bad" XI. We see the very 
significant effect that degradation of MWPC resolution in a single counter 
has on mass resolution. 

There are a variety of other measurement errors which contribute 
to mass resolution that we have included in our Monte-Carlo simulation, 
which have been described in detail by Wiedenbeck (1978). They include 
photoelectron (pe) statistics, electronic noise, and uncertainties in 
in-flight calibratious. The error in AE due to photoelectron statistics is 

= ■•'I ^ ■ . A tjrpical AE in D5+D6 for *®C is '**’ 1 GeV. With 

VlSOMev-^AE 

x/R« 0.6, Equation 3.B gives Cpe** 0.08 amu. To estimate the error in 
mass due to the finite statistics for calibrations we quadrature the 
0.15%, 0.25%, and 0.4% uncertainties from the Bevalac spatial response 
map, flight residual map, and PMT temperature correction, respectively 
(Chapter 2), to get « 0.5% uncertainty in both E5 and Ee. With AE = E5+Ee 

and EsW 2E5, we have )»uti«Uc3 “ 0.37%. Equation 3.8 
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then givea, for'- x/RsO.6, an^gtatiaUei'* Sach of these rough 

calculations is consistent with the more precise Monte-Carlo results. 


Three other effects which contribute to the mass resolution have 
not been included in the Monte— Carlo simulation; their magnitude has 
been evaluated instead by separate calculations and added in 
quadrature to the Monte-Carlo results to give the total expected mass 
resolution. The first of these is the possible error due to the assumption 
that E/L is linear in L. Its magnitude was estimated in §3.3 to be 0.05, 
0.07, and 0.11 amu for boron, carbon, and nitrogen, respectively. Second 
is the error we make in treating scintillator thickness variations as a 

correction to AL independent of x/R (§2.2.3). If is the percent 

w 


correction which should have been made to t, but instead was made to 
AE, then the mass we calculate is in error by 


»5M [AE 


3M 

dAE 


. 3M , (5t 
^ 3t ■' t 


I 




M 6tr ^ 

X 


- IJ. 


It is because the term in brackets varies with x/R that our treatment of 
thickness variations is not precise. That term, however, has an rms 

variation of only ~ 0. 11 for 0.4$ $ 0.9. From Table 2.1 ~ is « 0.6% for 

R t 


nitrogen, and slightly smaller for boron and carbon (as.suming thickness 
variations in D5 and D6 are uncorreiated). Thus our thickness error is 
<: 0.01 amu, and is therefore negligible. Finally, if the assumption that 
there is no spatial variation in the output of detector D7 beyond that 
meastired at the Bevalac (§2.2.3) is false, then part of the residual spatial 


% 
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variation assigned to the output of detectors D5 and D6 should in fact 
have been assigned to 07. The error in mass which results has two 
components whose effects add [this is because the calibration function 


fQ(sec 19 , 07 ) penetrating oxygen has 


afe 

a07 


« 1 . 


whereas the Tiass 


tracks (Figure 3.1) have 4^ < 0]. 

OL7 


The uncorrected mass for a’ ^vent at 


a particular position is 


aM 


aM 


M = Mo + rsE'-||i -(l-r)»vC55'iE-g^j, 


where Mq= actual mass of particle, r= fraction of the measured residual 
variation due to 07, and s= magnitude (%) of the residual variation (1.9% 
rms. Table 2.1). The factor under the square root comes from the 
assumption that 2 Eg and uncorrelated variations in signals Dg and 
Dfl (this factor is unity for nitrogen where Dg is not used). The (false) 

K dM 

correction we apply is vtD.55 ihai the error is 

™ = r=(E'||+vE:5FAE^) 

R R R 

We estimate r from Table 2.1 by assuming that the magnitude of the 

residual variation is proportional to that of the Bevalac variation, so 
1 _j. 57 

that = -tV» or r w 0.42. For boron and carbon we thus have 

r 4.2 

(averaged over 0.55^ ^ 0.90) 6M w 0.16 and 0. 19amu, respectively. For 

nitrogen, <5M w 0.28 amu (vb.55 -♦ 1, so that dM is independent of x/R), 


From Table 3.2 we see that the achieved mass resolution for all 
elements is bracketed by the calculation in which we assume nominal 


\ 
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MWPC XI resolution and 'that in which we assume degraded MWPC XI 
resolution. From the table it is also evident that there is a fundamental 
limitation to mass resolution because of multiple Coulomb scattering 
and Landau fluctuations of 0.15 to 0.22 amu which would dominate the 
total mass resolution if the MWPC resolution were the nominal ± 1 mm 
and there were no uncertainties in residual spatial variation. Improved 
MWPC resolution would obviously be desirable for all charges, as would 
an Independent calibration of passible residual spatial variations in 
detector D7. Other improvements in mass resolution could be obtained 
by using a detector material with a smaller 2^ than Csl to decrease the 
Coulomb scattering (Equation 3.8, although this could adversely affect 
other features: for example, there would be more nuclear interactions in 
the detectors), and improved statistics for in-flight calibrations. 
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Chapter 4 

Interpretetione of Meaenremente 

Having measured the isotopic compositions of boron, carbon, and 
nitrogen, we next interpret the measurements in Light of those 
processes — nucleosynthesis, galactic propagation, and solar 

modulation — by which they are determined. Before doing this, however, 
we must correct the measurements for the contamination introduced by 
the 5gcm"^ of atmosphere under which our experiment flew, as well as 
possible contamination due to the presence of the ’^10 g cm“^ of Csl 
(D0-D4) above our mass measurements at D5 (or D6 for nitrogen). 

4.1 Detector Interactions 

Our charge-consistency requirement (§3.2) eliminates most 
nuclear interactions in the Csl itself. It is possible, however, that 
neutron stripping reactions take place in the Csl above D5, altering the 
isotopic composition of a given element from its value at the top of the 
instrument. Based on the work of Lindstrom 9t al, (1975) we estimate 
that about 4% of the total inelastic cross section for a given particle in 
some mediixm can be ascribed to single neutron stripping. For the total 
inelastic cross section of a particle of mass number Ai in a medium of 
mass number A we use 

J. JL 

CT| = 52.20 mb [A ® + Al ® — b (4.1) 

where the overlap parameter is b = 1.189 exp [-0.05446 min (A,Aj)] 
[Hagen ei al. 19Y7; note, however, that b in this reference is given 
incorrectly as one minus the correct b used here (Ormes 1981)]. With 
Aw 130 for Csl we thus calculate, for example, that < 0,5% of will be 
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traxiffonned into in 10 g cm"^. We have therefore made no 
correction to the measured Isotopic compositions for this effect 

Because elemental abundances at the top of the instrument will be 
required for the atmospheric correction, we first use Figure 3.5c as a 
measure of the elemental abundances of particles which have undergone 
no charge-changing nuclear interaction before stopping in D7. The 
abundance at the top of the instrument for a given species is larger by a 
factor exp (0.96 norjx), where n» number density of nuclei in Csl and 
xsmean total thickness of Csl traversed for particles which stop in 
D7 n 52mm (the factor 0.96 is included assuming 4% of tX[ is due to 
neutron stripping). We list in Table 4.1 the measured elemental 
abundances from Figure 3.5c, the correction factors exp(0 96n(rjx), and 
the corrected elemental abundances at the top of the instrument. The 
fourth column gives the elemental abundances corrected for 
atmospheric secondaries and energy interval differences (sec following 
section), normalized to carbon. A comparison with another recent 
balloon measurement in the fifth column (Lund ut aL 1975) shows good 
agreement. 

4.2 Atmospheric Correction 

Based on independent measurements by our own 
pressure-transducers and those of the launch support group, the mean 
atmospheric depth at which the Instrument flew was 4.7± 0.2 g cm"^; the 
rms temporal variation in depth over the period during which data were 
taken was ~0.6 g cm"^. Because the atmospheric correction to our 
measurements tui is out to be quite small, temporal variations in depth 
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as well as the uncertainty In the mean depth have been neglected, and 
vre have simply treated the atmosphere as a 5 gcm'^ thick slab (includes 
a mean secant of 1.07) through which all particles must travel before 
encountering the instrument. 

The abundances at a given atmospheric depth and energy per 
nvciecn interval differ from those at the top of the atmosphere because 
of nuclear fragmentation and ionization energy loss. If fiCx.cn.cig) is the 
flux at depth x (g cm"®) of species i [atomic number Zi, mass Mi (pmu)] 
with energy per nucleon (MeV pmu"*) between Cu and Cj 2 , then the flux at 
depth x+6x is 

f,(x+dx,eii,ei2) * exp(=nc?i<5x)fi(x.c',i.e'i2) 

+ S [1 - exp (-noutix) ] fk(x,eiki.cik 2 ) (4.2) 

k>i 

where n » number of target (atmospheric) nuclei per gram of 
atmosphere, ai= total inelastic cross section of speck's i in the 
atmosphere (Equation 4.1), and tru,* partial cross section fcT species k 
to produce species i in the atmosphere. We have made the 
approximation that dx is sufficiently small that one can neglect tertiary 
and higher order fragmentations. The energy intervals are related by 

IW(e'ii) = R.i,fe,) + ^Sx. and = R.ir(=ii) + (-^ +^)-f 

similar equations relating e'i 2 and eu ,2 to ei 2 - H|^(c) is the range in air 

(g cm"®) of a proton of kinetic energy e. We have assumed that 

Sx 

fragmentations of species k occur at x + ~ . 

C 

We have developed a computer program based on Equation 4.2 
which take.*;, as input, abundances and energy spectra of nuclear species 




go 


at the top of the atmosphere and propagates these abundances in 
1 g cm~^ steps into the atmosphere. For the partial cross sections 7{]( we 
hnve used the semi-empirical proton-nucleus cross sections and scaling 
xiucAeus-nucleus cross sections of Silberberg and Tsso (1973a, 1973b, 
i977a, 19 l'7b, 1977c). (The energy dependences of the nucleus-nucleus 
cross sections are not well known: wc have used the asymptotic high 
energy values which are applicable down to at least 600 MeV pmu“^) We 
have used the tabulation by Barkas and Berger (1964) for the function 
As input to the program we have used the satellite measurements 
of Garcia-Munoz and Simpson (1979) for the elemental abundances of 
species with Z> B, and the calculation by Wiedenbeck (1970) lor the 
corresponding isotopic abundances. We have assumed that all species 
have the same energy spectra which we take to be that measured by 
Garcia-Munoz et al, (1977b) corrected to the level of solar modulation 
appropriate to August 1976 (Appendix C; specifically, we have assumed 
that df/dt = kc* with a = 0.57 for c< 220MeVpmu“' and a = -0.46 for 
c> 220MeVpmu"^). We have adjusted the elemental abundances of B, C, 
N, and 0 and the isotopic abundances of 6, C, and N (we have assumed 
that our observed 0 is all ^®0) at the top of the atmosphere until the 
abundances at 5gcm”^ are consistent with our measured elemental 
abundances (third column, Table 4.1) and isotopic composition 
(Table 3.1). 

We show in Table 4.2 the fraction of the abundance of a given 
nuclide at 5 g cm"^ atmospheric depth due to the various contributing 
species at the top of the atmosphere. For example, 07.6% of at 
5 g cm“^ comes from at the top of the atmosphere which undergo no 


Table 4.H Contributions to the Observed Abundances (Percent) 
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nuclsar interaction, 1.6% of results from nuclei which are 
stripped of a single neutron In 5 g cm**^ of atmosphere, and so on. That 
no more than 15% of the observed abundance of any of the nuclides in 
Table A.2 is due to fragmentation of heavier species suggests that our 
assumptions about the isotopic composition of species with Z 2 8 and the 
energy spectra of all species are relatively unimportant in making the 
atmospheric correction. 

We summarize in Table 4.3 the results of the atmospheric 
correction. The first column indicates the energy interval at the top of 
the instrument over which HEIST is sensitive to the given nuclide. The 
lower limit corresponds to the requirement that L 7 > 1 1 nlu (see 
Figure 3.2b), the Vipper limit results from the restriction that the 
vertical projection of the range in detector D7 must be less than the 
thickness of D7 (17 mm). The second column, from Table 3.1, gives the 
measured isotopic abundances. The third column indicates the amount 
by which an abundance at the top of the instrument must be increased 
to obtain its abundance at the top of the atmosphere. The fourth 
column gives the energy interval at the top of the atmosphere, and 
finally, in the fifth column, we tabulate the isotopic abundances 
corrected for both atmospheric interactions and the fact that the 
energy per nucleon intervals for isotopes of a given element are of 
slightly diflferent width. The overlap of the energy intervals is such that 
the fractional isotopic abundances vary by only ~0.01 as the assumed 
form of the spectrum varies from to in Table 4.3 we have 

simply assumed that the number of events is proportional to the width 
of the energy window, that is, a flat energy spectrum. 
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4.3 Comparisons with Other Measurements 

We show in Figure 4.1 our corrected results for the isotopic 
compositions of boron, carbon, and nitrogen, together with those of 
other investigators. Tabulated in Table 4.4 are the mass resolutions 
achieved by this work and others. Our mass resolution at boron is the 
best of those tabulated. At carbon and nitrogen our resolution is 
comparable to or better than others*, except .r the exceptionally high 
resolution achieved by Wiedenbeck cf aL. (1”" ''V Our relatively large 
uncertainties in the determination of fractional isotopic abundances are 
due almost entirely to limited statistics. Because of the excellent 
agreement between predicted and achieved resolution (Table 3.2), we 
feel confident that systematic effects have been properly identified and 
accounted for. 

Also shown in Figure 4.1 are smooth curves which represent the 
predictions of a model whose characteristics we now describe, 

4.4 Theory 

The cosmic ray abundances observed at earth differ from those at 
the cosmic ray source for two reasons. One, the particles traverse, on 
the average, ~5-6 g cm”^ of interstellar matter before being observed in 
the solar system, and are thus subject to nuclear fragmentation and 
ionization energy loss. Thus observations of the elements Li, Be, and B, 
which are bypassed in the normal processes of nucleosynthesis and are 
presumed to be absent at the cosmic ray source, are explainable entirely 
in terms of the spallation of heavier species. Second, as the cosmic rays 
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Ficure 4.1 

Comparisons of measurements. C: this work. □ : Garcia-Munoz et 
0.1 1977a, Guzig 1980. X: Wiedenbeck ef al. 1979. A: Hagen ot al. 1977. 
V : Buffington et a/.. 1970. O : Webber and Kish 1979, Webber et al. 1979. 
The dashed curves correspond to a galactic propagation model with no 
solar modulation. The solid curves correspond to a level of modulation 
4 = 300 MeV amu"^ for particles with M/Z=2. 
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enter the solar cavity, they suffer energy loss as they interact with the 
small-scale irregularities of the interplanetary magnetic field which is 
carrried radially outward from the sun by the expanding solar wind. One 
must therefore account for the rigidity-dependent effects of this 
process before interpretation of observed abundances is possible. 


4.4.1 Galactic Propagation 


If ni(it,t,e)de is the number density of particles of type i at position 
X and time t with energy per nucleon between e and e+de, then the 
continuity equation which describes cosmic ray transport in the galaxy 
is (e.g. Reames 1974) 


=V-(-qVn,)+i[^),n,]-^+q,+ Ei 


(4.3). 


The first term on the right is the divergence of a diffusion current «iVni, 
which is due to spatial gradients in nt; is the diffusion coefiicient. The 
next term is a "divergence" in energy-space, and accounts for changes in 
ni due to ionization energy loss. We lump into the term -nj/rj all 
"catastrophic" losses, such as nuclear destruction and escape from a 
confinement volume; the mean time for such a loss is tj. Finally, there 
are two source terms, the production of particles at the cosmic ray 
source, and the contribution from heavier species j which fragment with 
mean time tij into species i. 


The usual approximations which go into simplifying Equation 4.3 
are i) steady-state equilibrium (3ni/8t=0), and ii) homogeneous 
distributions of the ui (Vni=0). It is convenient to consider pathlengths 
A (g cm~^) rather than times r [related by h-p^cr with pss density 
(g cm"®) of interstellar gas], and fluxes J rather than densities n 
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(jKn^c/47r). Equation 4.3 becomes 


(4.4) 


where Ji» flux of soecies i in local intersteUor space. Z|. charge, mass 
of species i, » mean free path for catastrophic losses of species i, 
Aij s mean free path for species J to fragment into species i, and 
S(c)s specific ionization (MeV cm^ g“*) of a proton of energy c in 
interstellar space. is proportional to the source production of species 
i per gram of interstellar matter. The solution to Equation 4.4 requires 
the specification of boundary conditions. In the currently popular 
'Leaky Box" model (Coswik e< al. 1967), the galaxy is envisioned as a 
confinement volume whose boundaries are encountered many times by 
the cosmic rays before they escape, implying an exponential distribution 
of pathlengths with mean A«. The solution to Equation 4.4 is then 
(Rcnmes 1974, Meneguzzi ei al. 1971) 

where c* is related to c through the integration variable x, J* ® \ ^ . 

The mean free path for losses is related to A* through ~ 

where Aid mean free path for nuclear destruction of species i. 

Equation 4.4 thus allows, for a given set of source abundances and 
energy spectra [embodied in the Qi(c)], the propagation of those 
abundances and spectra to give the fluxes Ji in local interstellar space, 
assuming that the Aij (which are related to the cross sections cij) are 
known. 
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4.4.2 Solar Modulation 

Fluxes In local interstellar space differ from those observed at the 
orbit of earth because of the influence of the expanding solar wind to 
which is tied on interplanetary magnetic field. The effects of convection 
and diffusion ore embodied in a spherically syn^metric Fokker-Planck 
equation (Fisk 1974) 

7 1; ■* 1?^ J It 

where n(r) s number density of particles per unit kinetic energy interval 
at distance r from the sun, V^s solar wind velocity w 400 km sec'^ 
diffusion coefficient, T« kinetic energy, and a»(T + 2Eo)/(T + Eo), with 
Eo being the rest energy the species of interest. The solution of 
Equation 4.8 predicts, for a given local interstellar spectrum and 
diffusion coefficient, the modulated spectrum at r s l AU. It is convenient 
to characterize the level of modulation by a parameter $, which 
represents the mean change in energy per nucleon from outside the 

T ^ V 

solar cavity to 1 AU (Gleeson and Axford 1968): $ w r— j dr— where 

□ M^au « 

Rw 50 AU is the distance beyond which solar modulation is negligible. 
The diffusion coefficient « is typically taken to be proportional to velocity 
times a function of the particle’s magnetic rigidity, Psy/SEo/eZ (Fisk 
1974, Lezniak and Webber 1971). A particularly simple form is ««/?P, in 
which case $ is independent of energy and proportional to the particle’s 
charge to mass ratio. For the period 1974-1976, Garcia-Munoz et 
al. (1977a) show that $ « 220 MeV amu"^ Wiedenbeck and Greiner 
(1960) take $ = 325 MeV amu"^ for the period August 1978 through 
August 1979 to account for an increase in solar activity. As our 
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measurttments were made in Auguat 1978, we eatimate 

♦ w 300 Me V am^.“‘. 

It ahould be emphaaized that t repreaenta the averape energy loaa. 
and doea not apply to individual particlea. Urch and Gleeaon (1973) have 
calculated the diatributiona of energlea outaide the aolar cavity from 
which came particlea with aome fixed energy at 1 AU. We ahow in 
Table 4.5, baaed on their work, the energy intervale in interatellar apace 
corresponding to the FWHM of the distribution, for our own 
measurements as well as the satellite measurements by the University of 
Chicago group (Din Figure 4.1, Garcia-Munoz «t al, 1977a, Guzig 1900), 
and the UC Berkeley group (X in Figure 4.1, Wiedenbeck «< al. 1979). 

4.4.3 Model Calculations and Discussion 

The models of galactic propagation and solar modulation 
described above have been used to predict the observed isotopic 
abundances of B, C, and N at 1 AU. (We are indebted to M. Wiedenbeck for 
supplying the computer-generated results of the calculation which we 
quote here. See Wiedenbeck and Greiner 1980, and references therein, 
for a further discussion of the model used.) The elemental composition 
of the source is taken to be that calculated by Silberberg et al. (1976). 
The isotopic composition is assumed to be that of the solar system 
(Cameron 1973), except in the case of the element Ne, for which recent 
measurements (Greiner et al. 1979, Mewaldt et al. 1960) indicate an 
abundance of ®^Ne enhanced over the solar system value. A source 
spectrum of the form (c + 400MeV amu”^)”^ ® (Garcia-Munoz et 
al. 1977a) is used, and partial cross sections are those of Silberberg and 
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Tsao (1973a, 1973b, 1977a, 1977b). The interstellar gas is presumed to 
contain one He atom for every ten H atoms; the H density is taken to be 
0.3 cm*’®. The mean escape pathlength is taken to be A, * 5.5 g cm”® 
(Ormes and Freier 1976). 

We show in Figure 4.1 (together with the experimental points) the 
results of the calculation. The solid curves correspond to a level of solar 
modulation i = 300 MeV amu"‘, the dashed cr*,rves correspond to no 
modulation (local interstellar space). Some qualitative features are 
readily apparent. The very slight energy dependence in each of the 
unmodulated curves is almost completely washed out by solar 
modulation, the first-order effect of which is to increase the fractional 
abundance of the heavier isotope at lAU. This is due to the fact that ^ 
is slightly larger for the lighter isotope, and local interstellar 
abundances decrease with energy. Second, quite apart from the 
predictions of the model, the data are consistent with a constant 
fractional isotopic abundance, independent of energy, for each of the 
elements. (Disagreements between individual measurements, which are 
apparently statistically significant, might in fact be due simply to an 
underestimate of the errors.) 

Ideally one would like to use plots such as this to place constraints 
on certain parameters of the model, although it is clear from the size of 
our error bars that stringent limits would not likely result from our 
measu. ements alone. Instead we discuss, in the context of the data as a 
whole, the degree to which interpretations of these data are limited, 
both because of measurement errors and model inadequacies. 
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Because the predictions of the model depend critically on the 
cross sections used in calculating the secondary component of the 

observed flux (the term 2 i in Equation 4.5), we show in Figure 4.2 

J>1 

the semi-empirical values of the most important cross sections as a 
function of energy. [By a cross section such as or (^®0+p“* *®N) we rea’ly 
mean the sum of the cross sections or (^®0+p-*^®N) and a (^®0+p-* *®0), 
because p* decays almost immediately to *®N.] Also shown are the 
energy intervals in interstellar space (Table 4.5) appropriate to the 
satellite measurements by the University of Chicago group 

(Garcia “Munoz si al. 1977a, Guzig I960), the UC Berkeley group 
(Wiedenbeck si al. 1979), and this work. The onset of a significant 
energy dependence to the cross sections below ~400MeV amu“^ means 
that solar modulation effects become more significant for the low energy 
measurements (as is seen in the shape of the dashed curves in 
Figure 4.1). More important, most of these cross sections have not been 
measured at low energies, and the form of their energy dependence can 
only be inferred from other measured cross sections, the nuclear 
physics of which is expected to be similar (we shall quantify this 
somewhat in the discussion of nitrogen below). With these limitations in 
mind, then, we discuss what one can conclude from the measurements. 

Because boron is produced entirely by the breakup of heavier 
species in the interstellar medium, the locations of the curves in 
Figure 4.1a are, to first order, dependent only on the ratio of the 
relevant cross sections for producting and For both species, 

oxygen and carbon are the major contributors, accounting for more 
than half of the observed boron. That *^B dominates is due in part to 
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Fiipire 4.2 

Semi-empirical cross sections (Silberberg and Tsao 1973a, 1973b, 1977a, 
1977b). 

a) Cross sections for on protons to produce the indicated products. 

b) Cross sections for *^C on protons to produce the indicated products. 

If the interaction can proceed by way of an intermediate short-lived 
radioactive nuclide which decays to the indicated product, the 
corresponding cross section is included. The horizontal bands 
correspond to the energy intervals in local interstellar space appropriate 
to measurements by the University of Chicago group (Garcia-Munoz et 
al. 1977a, Guzig 1900), the UC Berkeley group (Wiedenbeck et al. 1979), 
and this work. 
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the oontribuiioxi from the reaction ^®C(p,np)^*C; subsequently 

decays to We take Figure 4.1a as an indication that the ratios of 
the cross sections for producing and ^*B are adequately described by 
the senii“'empirical formulae, and show little energy dependence, at 
least above ~400 MeV amu"^ 

Turning next to carbon (Figure 4.1b), we see that a solar-like 
composition (^^C/C -sO.Oll) for the cosmic ray source is consistent vrith 
the measurements. If our model of galactic propagation is substantially 
correct, a precise determination of the source abundance is 

obscured by the large (> 00%) secondary contribution to the observed 
flux. The measurements do rule out a hot (T>10®°K) equilibrium 
CNO bi-cycle as a major contributor to the carbon abundances, as 
pointed out by Guzig (1980), as in that case one would expect a ratio 
*^C/C approaching 0.2. 

The nitrogen measurements (Figure 4.1c) are the most indicative 
of some disagreement with the predictions of the model, at least at 
satellite energies. In order to avoid requiring a source ratio of ^®N/N 
which is significantly greater than the solar system value of 0.0036, 
Guzig (1980) has used an alternate set of cross sections for the reactions 
i6o+p_, I5f^ l®0+p-* which are scaled from measured cross 

sections of and Using this alternate set as an 

indication of the sensitivity of the model predictions to uncertainties in 
cross sections, we find that the unmodulated value of *®N/N (dashed 
curve, Figure 4.1c) increases by '*'18%, 10%, and 4%, at interstellar 
energies (Table 4.5) appropriate to the measurements by the U. Chicago 
group, the UC Berkeley group, and this work, respectively. This implies 
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an increase in the predicted modulated Talue of by '*'0.07, 0.04, 

and 0.02. At satellite energies, then, the predictions of the model are 
uncertain at the same level as uncertainties in the measurements, even 
if one considers cross section errors alone. Finally, we find that a 
weighted average of the difference between measured and predicted 
'®N/N for these three measurements is 0.06 ± 0.04, where the error 
includes both cross section uncertainties and measurement 
uncertainties, the former as given above. We note that modest increases 
in solar modulation, and, to a lesser extent in the value of A^, would 
further close the gap between predicted and measured values of *®N/N, 
especially at low energies. We conclude that the source abundance of 
is marginally consistent with the solar system value. 

In conclusion we point out that, although at present there is no 
evidence that the isotopic compositions of boron, carbon, and nitrogen 
at the cosmic ray source differ substantially from those of the solar 
system, other recent measurements (e.g. Mewaldt et al, I960) have 
indicated enhancements in the neutron rich isotopes of neon and 
magnesium. Improvements in statistics, especially at balloon snergies, 
and in parameters of the model, most notably the low energy cross 
sections, may eventually show that and are similarly enhanced. 
Emphasis in the future might well be placed on high energy, high 
resolution, and long duration balloon measurements. 
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Chapter 5 
Summary 

We have * described a balloon-borne instrument capable of 
measuring masses of individual cosmic ray nuclei using 
position-sensitive multiwire proportional coiinters for measurement of 
particle trajectory, and Csl scintillators for measurement of particle 
energy, and have discussed the calibrations of that instrument using 
both data from an accelerator and the balloon flight itself. 

The saturation properties of Csl have been studied in the context 
of analyzing the data for mass; it is found that the diflerential 
scintillation efl!iciency of Csl decreases, for boron, carbon, and nitrogen 
nuclei, by '^1575 as the specific ionization of a particle increases from 
~50 to 200 MeV mm”*. A maximum likelihood technique has been used 
to analyze the mass distributions of boron, carbon, and nitrogen. The 
achieved rms mass resolution varies from ~0,3amu at boron to 
~0.5 amu at nitrogen, consistent with a detailed theoretical evaluation 
of the factors which contribute. 

The measurements have been corrected for nuclear interactions 
which occur in the Csl detector material, as well as for the 
contamination introduced by the ~5 g cm”^ of residual atmosphere. We 
obtain *°B/B = 0.3318;}!, *®C/C = O.OGlgioi. and *®N/N = 0.42lg;}? (the 
errors being primarily statistical), at energies of ~280, 300, and 
330 MeV amu”*, respectively. The results are in agreement with other 
measurements, both at lower (~ 100 MeV amu"*) and similar energies. 


We have described a model of galactic propagation and solar 


no 


modtilation with which we predict the abundemces observed near earth, 
given a solar-like isotopic composition at the cosmic ray source. We find 
no evidence for source abundances which differ substantially from those 
of the soiar system. 
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Appendix A 
IfWPC Baokfroand 

We show in Figures A.la-c the distribution in y position os 
measured by MWPC Y2 for penetrating events selected according to 
signal D7. For fast particles the output of detector D7 is proportional to 
(Equation 2.2) so that the three figures correspond approximately to 
different charge regions. Although the number of events within each 
sector agrees with the results of a trajectory simulation, the obviously 
anomalous tendency for the distributions to peak at sector centers is 
present to about the same degree in all MWPC's, and it is clear that the 
peaking becomes more pronounced for higher Z particles. 

The use of sectors and delay line readout in our MWPC’s [see 
Wiedenbeck (1978) for a comprehensive description] means that a given 
counter is unable to infer a two-particle event, and that the counter will 
respond to such an event by recording i) the sector of the particle with 
the larger ionization and ii) position within that sector represented by 
an average of the positions of the two particles within their respective 
sectors. This suggests that the peaking of the Y2 distribution toward 
sector centers is associated with two— particle events. In particular, 
because the feature is more pronounced for higher Z particles, we 
hypothesize that knock-on electrons produced in either the 0.69 g cm”^ 
aluminum shell which encloses our experiment, or the counter gas itself, 
are above threshold an appreciable fraction of the time. 

Because we have four (three) independent measurements of the 
y (x) position we can calculate Xx Xj (§2.2.1, Figure 2.3) according to 
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Fifure 

Histogram of y position as measured by MWPC Y2 for penetrating events. 
A crude charge selection has been made based on signal D7. 
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the airaight-'lizie ilia to th« x and y trajectoriea. Large x valuea indicate 
knock’-’on electron contamination, amall yaluea indicate high quality 
trajectory determination. Eventa for which Xy^ 3 mm have been 
re-analyzed under the aaaumption that one of the four y meaaurementa 
ia aignidcantly worae than the other three and ia therefore not uaed. 
Becauae of the relative geometry of XI, X3, and X4 (Figure 2.2c), Xz ^ 
eaaentially an indication of the extent to which X3 and X4 are conaiatent, 
and ia not uaeful in determining whether or not XI ia "bad". We are 
therefore unable to eliminate eventa for which XI ia contaminated by 
knock~ona. 

To determine what the degraded XI reaolution of auch eventa 
might be, we have uaed the y trajectory data »Ve show in Figure A.2 a 
histogram of the quantity < ys> -ya, that is, the deviation of the 
measured ya from that calculated according to the results of the 
straight-line fit to the y trajectory (only events for which Yl,2,4 were not 
contaminated are used in the dt; Y3 is always used). The central peak 
corresponds to "good" Y3 events, but there is a broad background with 
FWHM w 100 mm, or crft< 40 mm. Similar analyses of < yj> -yi for 1=1, 2, 4 
indicate that the resolution of a coordinate contaminated by knock-on 
electrons is a « 20-40 mm. 
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Figure A.2 

MWPC background as shown by a histogram of the quantity < ys> 
The vertical scale is logarithmic. 
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Appendix B 

Multiple Coulomb Scaiteriuf in Thick Oetecton 


To derive a form for the mean square multiple Coulomb scattering 
angle when the detector thickness x is large enough so that the particle 
velocity changes significantly in traversing it, we begin with Equation 3.6 
in differential form: 


da| X k 


Z^dx 


where k is Independent of velocity. With dx = and the fact that 

dE/dx 

scatterings in successive layers dx are independent processes so that the 
variances can be integrated, we have 

•t 

dx 

where E is the initial energy and E* the final energy. Since 

iS dE 


dx 


(/J,Z)w constant the quantity 

r 




dx 


Z« 


is to a good approximation 


O HF 

independent of both ^ and Z, so we treat it as the constant • 

In addition, the particle momentum p-y/JMmpC is related to the 
energy by (pc)^ = + 2 M mpC^E (mp is the proton rest mass). Thus 

k(mpc2)2 - = 


(t| = 




^ dE (£2 + aMmpC^E)-^ 


which integrates to 


X 


as: 


2M^2^(/3.1) 


In [■ 


2E 




(i+r)E' 7+1 


dx 


since E = ( 7 -l)MmpC®. 


IIB 


Appendix C 

Energy IfpectnrUsed in the Atmoapheric Correction 

Measurements by the Chicago group of the cosmic ray differential 
energy spectra (Garcia-Munoz at at. 1977b) indicate that the relative 
elemental abundances are, to a good approximation, independent of 
energy, at least over the energy interval 200-700 MeV amu"^. We have 
approximated their measurement of the oxygen spectrum by 

J(e)«k(-i)* (C.l). 

Co 

where Cq = 300 MeV amu■’^ and a = 0.35 for c< cq* ® * “0.77 for e> cq. 
Because the level of solar modulation appropriate to their 
measurements was less than that for the HEIST flight (§4.4.2), we have 
performed the following calculation which predicts the spectrum for 
1978, 


The flrst-crder effect of the expanding solar wind on the cosmic 
rays is to reduce the value of their energy per nucleon outside the solar 
cavity to that at 1 AU by an amount $ . If Jis(c) is the differential energy 
flux outside the solar cavity in local interstellar space, then the flux at 
1 AU ('Torce-field" approximation, Fisk 1974) is Jiau(c) = s($,c) Jis(c+$) , 

where s($,c) = ^ Here Mq = 938 MeV is the proton rest 

(e+4'+Mp)2-M| 

energy. If we assume that the local interstellar flux does not change 
with time, then we can relate the fluxes at 1 AU, iJiau 2^1 AU» 

characterized by the levels of modulation and respectively, by 

.J.Au(t) = if*- *«('+♦ 1~*2) • 


With 2 J 1 AU equal to the measurement by Chicago given above. 
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220 McV amu"^ and 300 MeV amu“‘, we find that iJiau t>e 

approximated by a form similar to Equation C.l, except that 
Co * 220 MeV amu“^, and a » 0.57 for e< Cq, a = -0.46 for c> c©. This, 
then, is the form we use in making the atmospheric correction. 
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